Strain induced correlation gaps in carbon nanotubes 
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We calculate the change in the correlation gap of armchair carbon nanotubes with uniaxial elastic 
strain. We predict that such a stretching will enlarge the correlation gap for all carbon nanotubes 
by a change that could be as large as several meV per percent of applied strain, in contrast with 
pure band structure calculations where no change for armchair carbon nanotubes is predicted. 
The correlation effects are considered within a self-consistent Hartree-Fock approximation to the 
Hubbard model with on-site repulsion only. 

PACS numbers: 62.25.+g, 71.10.Pm, 71.20.Tx 
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Carbon nanotubes (CNT) have many extraordinary 
electronic and mechanical properties In particular, 
band theory predicts that CNT are either metallic or 
semiconducting depending on chirality, i.e. in which 
direction a graphite monolayer is "rolled up" into a cylin- 
der forming the tube |2J. Probing CNT with a scanning 
tunneling microscope, this metallic or semiconducting 
behaviour could be tested experimentally 0- GJ- In 
addition, CNT can sustain large mechanical strains and 
can be deformed elastically up to bendings of order 19° 
which corresponds to a strain along the tube of 5.5% 
Experiments and numerical calculations indicate a large 
Young modulus of order lTPa ||J. The interaction of 
mechanical and electronic properties has been studied at 
room temperature in two experiments 0, Q where it has 
been shown that uniaxial stress can change dramatically 
the electronic structure of CNT. In these experiments 
CNT have been suspended between two metal electrodes 
on SiC>2/Si substrates. To apply uniaxial stress along 
the tube the tip of an atomic force microscope (AFM) 
was used. The tip was lowered to push at the center of 
the CNT. The AFM allows to measure simultaneously 
the deflection and the conductance of the CNT. The 
strain can be defined as a — [V4<5 2 + I 2 - /] // where I 
is the suspended length of the tube and 5 is its vertical 
deflection. In the first experiment |f| it was shown that 
changing the strain from 0% to 3.2% let the conductance 
of a metallic CNT drop by two orders of magnitude. 
Both the mechanical and the electronic properties were 
observed to be completely reversible. More recently 
it was demonstrated Q that not only metallic CNT 
become less conducting when applying stress but also 
inversely that some samples modified their behaviour 
from semiconducting to metallic. These experiments 
show convincingly that uniaxial stress applied to CNT 
changes their electronic properties. 

Theoretically, the effect of strain on the electronic 
properties of CNT has been studied in band structure 



calculations, cither analytically, using a tight -binding 
approach, or numerically by density functional theory 
[E IS EES IH El- In particular it has been shown 
that, depending on chirality, uniaxial stress can increase, 
decrease or not alter the band gap. 

In this paper we include the effect of electron-electron 
correlations in these calculations. We compute in 
detail the gap as a function of applied strain and we 
compare our results to the one-electron band structure 
calculations in the literature. The calculations are 
carried out within a Hubbard tight -binding model using 
the self-consistent Hartree-Fock (H-F) approximation. 
It has been argued by the authors ^3 that the charge 
gap of CNT at half-filling is well described within this 
approximation. 

In this work we consider single-walled CNT at half- 
filling. In a first approximation, single walled CNT can 
be thought to be a rectangular graphite monolayer with 
the appropriate boundary conditions (fig. They can 
be classified by their chirality vector Ch — na\ + maj, 
where ai and are the basis vector of the honeycomb 
lattice, while n and m are integers with m<n||. Ch de- 
termines into which direction the graphene layer is rolled 
up. We model CNT with a Hubbard model with nearest 
neighbour hopping between tt orbitals only, and on-site 
interaction at half-filling: 
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a is the spin index and i , j sum over the sites of a 
rectangular honeycomb lattice with periodic boundary 
conditions. cj CT (ci a ) are the fermion creation (annihila- 



tion) operators and rii a 



The hopping integrals 
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tij are restricted to nearest neighbors and in general 
they can have different values in every hopping direction, 
say t±, tiu and t 2 \\ (cf fig. P). These three hopping 
amplitudes and the one-site interaction strength U are 
the four parameters entering into the model. Choosing 
some values for those parameters we can determine the 
charge gap through a H-F calculation for a given (n, m) 
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FIG. 1: A rectangular honeycomb lattice of a (4, 2) CNT, i.e. 
Ch = 4ai + 2d2- The basis vectors are chosen to be a± = 
o/2(V3, 1) and a 2 = a/2(V3, -1). a = 2.49A is the lattice 
constant for CNT. The chirality vector Ch and the vector T 
define the quasi-one-dimensional unit cell. Ch describes the 
circumference and T is oriented parallel to the tube. 



CNT. All the details of the H-F approximation can be 
found in reference |14| where also its validity has been 
discussed. 

If there is no applied strain we set in our calculation 
all the hopping amplitudes equal to a certain value 
to- This is of course only an approximation. Forming 
a CNT from a graphite sheet, where this is true, we 
change the hopping integrals due to effects of curvature. 
It has been shown 0, ^| that there is a curvature 
induced gap due to a-ir hybridization for all CNT 
except armchair CNT. This gap is at least of order 10 
meV. Thus armchair CNT are the only CNT where 
band structure calculations predict a metallic gapless 
behaviour. They correspond to chiralities with n = in 
or equivalently to a chiral angle 6 = 7r/6 (cf fig. QJ. For 
armchair CNT the curvature induces only two different 
hopping amplitudes, namely t± and t|| = t\\\ = t 2 ||. Due 
to this symmetry, no gap is opened from the band point 
of view. However we will show that there will be a gap 
induced from electron-electron correlation effects for all 
CNT, even if they are of armchair type. We concentrate 
our discussion on armchair CNT. The same effects are 
apparent in CNT of other chiralities but in that case 
these effects are less visible as there is already a gap in 
the one-electron band structure. 

Following references 0, Q we use Harrison's phenomeno- 
logical law to relate the hopping parameter to of the 
undeformed CNT to the ones of the elastically deformed 
CNT 16] U = to (rl/ri) 2 where r' t and respectively, 
are the bond vectors before and after the deformation 
and i =_L, lii, 2n. Projecting these vectors along the 
directions of T and Ch, we can write for an elastic 
uniaxial strain along the tube: r,*T = (1 + cr)r' iT and 



r iC h = (1 — VG ) r 'ic h wnere v is the Poisson ratio, to can 
be estimated from ab initio calculations to be 2.4 eV [lj 
and the Poisson ratio has been computed numerically [lj 
and measured experimentally for graphite to be about 
v = 0.2 ^^|. It is difficult to get an estimate for the 
on-site repulsion U of atomic carbon. In the literature 
values between 5 and 12 eV are suggested [l9l I20I l2lj . 

To study now the correlation effects in the situation 
described above, i.e. stretching the tube by application 
of uniaxial strain, we rely as in previous work [l4j on the 
H-F approximation. Based on the observation that the 
H-F calculations reproduce correctly the functional de- 
pendence of the energy gap in the one-dimensional Hub- 
bard model at half filling, we believe that it gives also 
reliable results for the ga p of CNT. 

We have shown in Ref. [l4[ that the H-F calculations 
for the charge gap of armchair CNT give the following 
scaling law 

E g /t\\ = c/n exp {-an{t\\/U - t^/U?/)} (2) 

with c = 1.01, a = 5.44 and U^ F = 2.23t a , the critical 
value to open a gap in the two-dimensional honeycomb 
lattice. Since graphite is a semi-metal U is expected 
to be smaller, but close to this value, i.e. U < U^ r F ' . 
For such values a correlation gap of order 10 meV is 
present for CNT of small diameter. We found that the 
scaling law (0) is still valid when a uniaxial stress a is 
applied, which produces a change of the ratio in 
the hamiltonian from 1 to larger values, tx/tu and 
a are related by Harrison's formula in a non-linear way. 
The prefactor c depends only slightly on the change in 
t±/tn. The variation of a and U^ r F on the applied strain 
is much more significant and it is shown in figure [21 We 
observe that U^. F increases with applied strain. From 
equation (J2J we can see that this would imply that the 
charge gap diminishes since a U far from the critical 
value disfavours a large gap. However the simultaneous 
decrease of the parameter a overcomes this tendency 
and when both effects are taken into account, the charge 
gap increases approximately linearly with the strain 
with a slope that depends on U. Note that U^ F > 
indicates that there is a metal-insulator transition at a 
finite value U — U^ F ■ This is due to the fact that the 
density of states in the honeycomb lattice is vanishing 
at the two Fermi points as p(e) cx |e|. In the limit 
t±/t\\ — > we get the one-dimensional behaviour where 
any infinitesimal electron-electron interaction can open 
a gap. The other limiting point is fj_/f|| = 2 where 
the two Fermi points collapse into one and a band gap 
appears for values t±/tu > 2 which makes the system 
insulating already for U = 0. 

To compare our results to experiments we plot the 
variation of the gap for armchair CNT as a function of its 
size when a strain of 1% is applied (cf figure 0). The on- 
site interaction U is set to 2to- As the hopping amplitude 
to has a large value of 2.4 eV neither the correlation gap 
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FIG. 2: Numerical results for the dependence of the parame- 
ters a and U^ r F in equation on t±/tn: (a) a is extracted 
from fitting the H-F results to the scaling law J2J . (b) U^. F is 
calculated by evaluating xo = J e<0 de^r, the bare staggered 
static susceptibility, and where p(e) is the tight-binding den- 
sity of states of the honeycomb lattice. In H-F theory we 
have Uc r F = Xo 1 ■ 



FIG. 3: The variation of the correlation gap for different 
(n, n) armchair CNT is plotted when a strain of 1% is ap- 
plied (U — 2to). The upper plot shows the variation of the 
gap E g (a — 1%) — E g (a — 0) as a function of n. n is pro- 
portional to the tube diameter. The lower plot shows the 
same variation as a function of the original correlation gap at 



zero applied strain E g o = E g (a = 0). 
different armchair CNT, n = 5, . . . , 26. 



E g o is calculated for 



at a = of order 10 meV nor its variation of order meV 
per percent of applied strain can be neglected. As an 
example we look at the values for a (10, 10) CNT. We can 
read off from figure [3] that the energy gap at zero strain 
is 14 meV which corresponds to a temperature of 160 K 
and that a strain of 1% induces a rise of 4 meV in the gap. 

How does this compare to the effect of strain to the 
band structure? For small strains the following formula 
has been derived from a tight -binding calculation [8j: 

dE 

-j-^ =sgn(2p+ l)3to (1 + f) cos 39 (3) 
da 

where p G { — 1,0,1} is defined by the equation 
n — m — 3q + p (q is integer) . This formula has been 
used to interpret the experimental results in Ref. |7j. It 
follows from it that the change of the band gap with ap- 
plied stress can be either positive or negative, depending 
on the value of q, or in other words the chirality. The 
maximum variation is achieved for zig-zag CNT (chiral 
angle zero) and is about ±85 meV/%. The maximal 
variation of 85 meV/% is one or even two orders of 
magnitude larger then the variation of 4 meV/% which 
we derived from electron-electron correlation. However, 



for armchair CNT (8 = ir/6) equation © and ab initio 
calculations predict that no gap opens up with applied 
strain. Then correlation effects are the only reason 
why one should have a gap and this gap increases by 
applying uniaxial strain at a rate of several meV/%. In 
our previous example of the (10, 10) CNT, we have seen 
that one percent of strain can change the correlation 
gap by about 30% of its original value. 

To summarize, for semiconducting CNT with large 
band gaps the electronic structure at half filling is well 
described within band theory and correlation effects give 
only small corrections. But for CNT with a small gap 
(of order meV) correlation effects cannot be neglected. 
This is especially true for armchair CNT where no band 
gap at all is predicted but they should develop a mea- 
surable gap, induced from correlations alone, if sufficient 
pressure is applied. This conclusion is illustrated in fig. 
01 We plot the strain which is necessary to open a gap 
of Eg = 50mcV as a function of the on-site repulsion 
U for a (10, 10) armchair CNT. We see that if U is not 
too far from U cr this quite large gap would be realisable 
experimentally and could be seen in low temperature 
experiments. 
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FIG. 4: Minimal strain a min necessary to open a gap of 
50meV as a function of the on-site repulsion U. The calcu- 
lation was done for a (10, 10) armchair CNT. 
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